APS/123-QED 



Antiferromagnetic superconductors with effective mass anisotropy 

in magnetic fields 
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We derive critical field Hc2 equations for antiferromagnetic s-wave, d^2_y2-w&ve, and dxy-^scve 
superconductors with effective mass anisotropy in three dimensions, where we take into account 

(i) the Jaccarino-Peter mechanism of magnetic-field-induced superconductivity (FISC) at high fields, 

(ii) an extended Jaccarino-Peter mechanism that reduces the Pauli paramagnetic pair-breaking effect 
at low fields where superconductivity and an antiferromagnetic long-range order with a canted spin 
structure coexist, and (iii) the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO or LOFF) state. As an 
example, experimental phase diagrams observed in organic superconductor K-(BETS)2FeBr4 are 
theoretically reproduced. In particular, the upper critical field of low-field superconductivity is well 
reproduced without any additional fitting parameter other than those determined from the critical 
field curves of the FISC at high fields. Therefore, the extended Jaccarino-Peter mechanism seems 
to occur actually in the present compound. It is predicted that the FFLO state does not occur 
in the FISC at high fields in contrast to the compound A-(BETS)2FeCl4, but it may occur in low- 
field superconductivity for s-wave and ci^2_y2-wave pairings. We also briefiy discuss a possibility of 
compounds that exhibit unconventional behaviors of upper critical fields. 
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I. INTRODUCTION 

Recently, magnetic-field-induced superconductivity 
(FISC) has been observed in organic superconductors 
A-(BETS)2FeCl4 and K-(BETS)2FeBr4 [ES, 3, where 
BETS is bis(ethylenedithio)tetraselenafulvalene. In these 
salts, localized spins on Fe'^+ exhibit antiferromagnetic 
long-range order at ambient pressure at low temper- 
atures. The FISC in these compounds is considered 
to be due to the Jaccarino-Peter mechanism d, Q, 
where the localized spins are aligned uniformly at high 
fields. Konoike et al. have observed in the compound 
K-(BETS)2FeBr4 that superconductivity coexists with 
the antiferromagnetic long-range order in a low-field re- 
gion around the zero field [s']. They have fitted the ex- 
perimental phase diagrams by Fisher's theory Jj] based 
on the Jaccarino-Peter mechanism. The resultant up- 
per critical field of low-field superconductivity is much 
smaller than their experimental data. They have sug- 
gested that the reason for the discrepancy is that the 
Jaccarino-Peter-Fisher theory does not take into account 
the antiferromagnetic long-range order at low fields 

In recent works, one of the authors has extended 
the Jaccarino-Peter mechanism to antiferromagnetic su- 
perconductors with canted spin structures in magnetic 
fields 0, The canted spin structure generates the 
ferromagnetic moments that create exchange fields act- 
ing on the conduction electrons through Kondo interac- 
tions. If the Kondo interactions are antiferromagnetic, 
the exchange fields partly cancel the Zeeman energy. As 
a result, the Pauli paramagnetic pair-breaking effect can 
be largely reduced, and the upper critical field can ex- 



ceed the Pauli paramagnetic limit (Chandrasekhar and 
Clongston limit) Q. This mechanism occurs even in 
the presence of the orbital pair-breaking effect _&]. We 
call this mechanism an extended Jaccarino-Peter mech- 
anism in this paper. Since the canted antiferromag- 
netic phase occurs in the compound K-(BETS)2FeBr4 for 
H II c [l,[l3|j we apply the mechanism to this compound. 

In the compound K-(BETS)2FeBr4^he FISC has been 
observed both for H || c and H || a The phase dia- 
grams for H II c and H || a are rather different, and it is 
attributed to the anisotropy of the Fermi surface and the 
Kondo interactions between the localized spins and the 
conduction electrons. We take into account the Fermi 
surface anisotropy by effective masses. 

The effective-mass model was introduced in Ginzburg- 
Landau equations by Ginzburg [ll|, [l^ . Hohcnberg and 
Werthamer [31 pointed out that detailed structures of 
the Fermi surface affect the upper critical field. Rieck 
and Scharnberg [l3| and Langmann [l^ obtained gen- 
eral equations for arbitrary Fermi surfaces. Xu et al. [l6[ 
and Kim et al. [13 calculated the upper critical fields of 
mixed d-wave and s-wave superconductors with effective- 
mass anisotropy. Recently, Kita and Aral ^Ts'l have for- 
mulated an equation for the upper critical field, taking 
into account the Fermi surface anisotropy and the gap 
anisotropy on the basis of the Rieck and Scharnberg the- 
ory [14]. They have performed the quantitative calcula- 
tions of the upper critical fields for type-II superconduc- 
tors Nb, NbSe2, and MgB2 using Fermi surfaces obtained 
by first-principles calculations [l8[ . 

A theory of the upper critical field for layered super- 
conductors has been proposed by Lebed and Yamaji 
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and developed by Lebed and Hayashi (20j. They have 
found that when the layer spacing is large the upper crit- 
ical field exhibits a reentrant transition or an enhance- 
ment at low temperatures in the quantum region, due to 
an effect of dimensional crossover induced by the mag- 
netic field [T9|. In the compounds K-(BETS)2FeBr4, how- 
ever, since the upper critical field of the low-field super- 
conductivity did not exhibit either a reentrant transition 
or an enhancement in the experimental phase diagrams, 
the dimensional crossover does not seem to take place. 
Therefore, from a phenomenological consideration, we 
use the effective-mass model as an approximation instead 
of strict equations in Refs. and [23] ■ The effective- 
mass model is adequate in the Ginzburg-Landau region 
for layered superconductors. 

In this paper, first, we derive critical field equations 
for s-wave, dx2_y2--wave, and d^y-^ccve superconductors 
with effective-mass anisotropy in three directions, tak- 
ing into account both orbital and paramagnetic pair- 
breaking effects. Secondly, we take into account the ex- 
tended Jaccarino-Peter mechanism. Lastly, we reproduce 
the phase diagrams of K-(BETS)2FeBr4 including both 
the FISC and low-field superconductivity. 

We also examine the possibility of the FFLO state. 
The FFLO state has extensively been studied [2l| since 
pioneering works by Fulde and Ferrell, and Larkin and 
Ovchinnikov [2^. The state is taken into account by 
an extension of the BCS mean-field theory to include 
the possibility of finite center-of-mass momenta q's. In 
this study, we adopt a model in which q y H is as- 
sumed following Gruenberg and Gunther [23|, since we 
consider the situation in which substantial orbital pair- 
breaking effect is present. In the organic compounds 
A-(BETS)2FeCl4, the possibility of the FFLO state in 
the FISC has been pointed out by Uji et al. [H and Bal- 
icas et al. [1| , and also examined theoretically 0, H^l ■ 
The shape of the phase boundary of the FISC is well 
reproduced by taking into account the FFLO state [25j . 
Tanatar et al. have also argued that the FFLO state 
may occur in A-(BETS)2GaCl4 from their experimental 
data [26| . Recently, the quasi- two-dimensional heavy- 
fermion superconductor CeCoIns has been believed to 
exhibit the FFLO state '281, IM H IMI . Adachi and 
Ikeda [S^ and Won et al. [33] have calculated the criti- 
cal fields for dj.2_y2-waNe pairing taking into account the 
FFLO state in connection with CeCoIns. 

This paper is constructed as follows. In the next sec- 
tion, we extend the theory of the critical field to the 
systems with anisotropic Fermi surfaces. The s-wave, 
(ia.2_„2-wave, and d^^jj-wave pairing superconductors are 



examined. In Sec. IIIIl we take into account the extended 
Jaccarino-Peter mechanism 0, g. In Sec. EVl we ap- 
ply the present theory to the organic superconductor 
K-(BETS)2FeBr4. We compare experimental and theo- 
retical phase diagrams. The final section is devoted to 
the summary and a discussion. 



II. UPPER CRITICAL FIELD IN SYSTEMS 
WITH EFFECTIVE-MASS ANISOTROPY 

In this section, we derive critical field equations for s- 
wave, d3;2_j,2-wave, and d^jy-wave superconductors with 
an effective-mass anisotropy. We consider the electron 
system with the energy dispersion 



e(p) 
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2mz 2miy 



(1) 



and the pairing interaction 
n' = j d^r J d3rV|(r)V;T(r)T^(r-r')^j(r')Vi(r') , (2) 

where we have defined p — (pi,P2,P3) = {Px,Py,Pz) and 
introduced anisotropic effective masses rrix — mi, my — 
7712, and ruz = rn^. In the magnetic field H = B = rot A, 
the Hamiltonian is written as 



Ti. — Ti.0 + Ti-m + 'H' 



(3) 



where 
Ho = 



E 



d'rV4(r) 

X E — 

^ 2777^ 
1^=1,2,3 

nm = - I d^r fi^n 



-ih- 



^?/'i(r)cr^^-V<T'(r) 



Here, il)a{'^) denotes the field operator which annihilates 
an electron of spin ct at a point r = (x, y, z) = {xi, X2, x^). 
We have defined the electronic magnetic moment /io = 
— 5cMb/2 with the Bohr magneton /ib = ft|e|/(2777c) and 
the g factor of the conduction electrons cjc- We use the 
units such that c — = h = \ unless it is explicitly 
expressed. 

We define the axis in the spin space along the di- 
rection of H. We should note that the axis does not 
necessarily coincide with the z axis of the electron coor- 
dinate depending on the direction of the H. Therefore, 
we have 



E 



d^Y cr/7V4(i-)^/'<T(r) 



(5) 



where h = fi^H with H — |H|, and a = +1 and —1 for 
up and down spin states, respectively. 

We apply the BCS mean-field approximation to the 
interaction Hamiltonian ([2]) as 



H' = - d^v d-^r 



Ax(r,r')7A|(r)7^|(r') 



+ A*(r,r')Vi(r')V'T(r) 



(6) 
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defining the order parameter 

A_.(r, r') EE -Vir - r')(V'-.(r')V'<x(r)) . (7) 

We define tlie center-of-mass coordinate R = (r + r')/2 
and the relative coordinate p = r — r', and redefine the 
gap function as A_o.(R, p). By following the procedure of 
Refs. t34i] and we obtain the linearized gap equation 

A-„(r, p) 

xF(p-p')GW(-p',_zc.„) (8) 

exp(ip -pjCL^lp ,-lUJn) 

X exp(ip • n)A_CT(r, p') 

near the second-order phase-transition point. We could 
examine both orbital and paramagnetic pair-breaking ef- 
fects with Eq. dH]). In the derivation of Eq. it is 
assumed that the spatial variation of vector potential 
A is sufficiently slow. Therefore, the order parameter 
becomes a slowly varying function of the center-of-mass 
coordinate. We have defined 



After simple algebra using Eq. ((TO)) . we obtain 



2e 

n = -ihV A(r) , 



(9) 



the Matsubara frequencies ojn — (2n-l- l)TrT with integer 
n = 0, ±1, ±2, . . . , and the free-electron Green's function 



G(°)(p. 



1 



(10) 



iuj„ — e(p) + crh + ^ 

We introduce p = {Px,Py,Pz) = {Pi,P2,P3) with 
Pi, — {m/m^y/'^p^ so that the dispersion relation be- 
comes isotropic as 



^9 



1 



in the p space, where m = {nixmymzY^^ . We also intro- 
duce f = {xi, X2, xs) = {x,y,z) with Xi, = {■m^/mY^'^Xi, 
and d/dxu = {jh/m^Y/'^d / dx^, so that p • r = p • f . We 
also define the operator 11 with 



n 



n 



(12) 



where p = (pi,p2,P3) and = {rriv/rhY/'^ p^- We will 
calculate the explicit form of 11 afterward. 

We replace |p| and |p'| in ^(p — p'), A_(j(r,p), and 
A_cr(r, p') with the Fermi momentum since electrons 
with momenta p and p' near the Fermi surface mainly 
contribute to the gap equation. Therefore, we write 
^(P-P') and A_cr(r,p) as V{p,f>') and A_^(f,p), re- 
spectively, with unit vectors p = p/|p| and p = p/|p|. 



A-.(f,p) 

\0Jn\<^Y) 
OO 



'0 

' t 



X exp 



2i 



sgn(tj„)vF(p ) ■ n 



A-a(f,p') 



(13) 



where we have defined vf(p) = p/m and the de nsity 
of states at the Fermi energy A^(0) = rhpF / {2Tr'^h'^) [3a |. 
Since Eq. (fT^ does not depend on the spin value a, we 
omit the spin suffix from now on. By redefining p' = 
tvF(p')/2 = tvpp'/2 and summing up the Matsubara 
frequencies, the gap equation (jl3p acquires the form 

A(f,p) = -M d'p'Vip,p')- 



X cos 



p' s\ah.{2'KTp' /vy) 
A(f,^')- 



^/.-fp'.n 



(14) 



It is obvious that if we rewrite the integral variable p 
as p, the linearized gap equation (fM]) is similar to those 
obtained by many authors [13, [ll Ijj, 124, l35|, l36|, isj. Ill] . 
The differences between Eq. (fT4|) and the equations ob- 
tained so far for isotropic systems are that the vec- 
tor potential A is scaled as Ai, — (my/m)^/^yly in n, 
and that the pairing interaction l^(p,p') is deformed in 
the p-space anisotropic superconductivity. The former 
scale transformation in A has been studied in Ginzburg- 
Landau theory [11, 12, 39] for s-wave pairing. For s-wave 
pairing, since T^(p,p') = V{p,p'), Eq. (fT4|) is exactly re- 
duced to the equations for the system with the isotropic 
Fermi surface, except the scale of A^^ which results in 
the scaling = {my/rfiY^'^Hy. Therefore, the only dif- 
ference due to the Fermi surface anisotropy is that the 
critical field in the v direction is enhanced or reduced by 
the factor [mlm^Y^'^ ■ The latter deformation of y(p, p') 
occurs because l^(p, p') ^ V^P, p') for anisotropic super- 
conductivity. We take into account the deformation of 
the pairing interaction in the p space. These changes af- 
fect the mass anisotropy dependence of the critical field. 

Now, we derive a more explicit form of the upper crit- 
ical field equation. We take the z axis in the direction 
where the effective mass is the largest. The compounds 
with layered structures with large layer spacing are ap- 
proximately described by the models with to^ 3> m^;, rriy. 
For example, in the application to the organic supercon- 
ductor K-(BETS)2FeBr4, we take the x and y axes along 
c and a axes of the compounds, respectively. We study 
a quasi-two-dimensional superconductor under magnetic 
fields parallel to the layers (H || a; or H || y) hereafter. 
We derive the critical field equations only for H || x, be- 
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cause those for H \\ y are obtained by exchanging nix 
and iriy. 

We assume the uniform magnetic field H in the —x 
direction paraUel to the layers, which is expressed by a 
vector potential A = {0,0, —Hy). Then, we obtain the 
explicit form of 11 as 



„ , 2|e| 

n = -^ri—, -in—, -ih— 

ox oy oz 



Hy 



from Eqs. ([9]) and (fT2|) . We define the differential opera- 
tors 



1 



4 



/2k. 
1 



/2k. 



in. 



in. 



(15) 



which satisfy the bosonic commutation relations, where 

m 2m ^ ^^^^ 



The factor mj ^myVfiz is due to the effective-mass 
anisotropy. 

We consider pairing interactions of the form 



^^(P,P') = -.9a7a(p)7a(p') 



(17) 



where 7q(p) denotes the symmetry function of symme- 
try a and a = s, dr^'i_yi, dxy and so on. The symme- 
try function 7q(p) is normalized so that ([7a(p)]^) = 1, 
where the average on the Fermi surface (• • • ) is defined 
by (•••) = ! dn^p{0,6,^){- ■ ■)/ S dnpp{0,e,^), where 
p(0, 0, tf) denotes the angle-dependent density of states 
on the Fermi surface. The gap function is proportional 
to the symmetry function 7q(p) and expanded by the 
Abrikosov functions 0n(f) of the Landau- level indexes 
71 = 0, 1, 2, ... , as 



A(f,^) = ^A^7„(p)0„(f)e^«^* 



(18) 



where we have defined 

2wf J p"' smh(27rJ p/uf) 



T 



\ _ g2pWD/5F 



2vp J p-^ sinii(2TT 1 p / v-p) 



^i{n' —n)(p 



J = l 



(22) 



and 



4^2,(p,0~)^exp 



[-^p^sin^fl- 



min(n.n') 

E 

*:=0 



p^ sin^ 



k+j 



\J n\n'\ 



kl{k + 2j)![min(n,n') - fc]! 



X cos 



2p,, , qvF . 

— (h-\ cos( 

vf 2 



(23) 



More explicit expressions of the matrix elements D^^,, 
D^^r"^ , and -D^^" are given below. Defining the zero- 
field transition temperature Tj"'' , the linearized gap equa- 
tion (PO)) is written as 



log 



T 



A" = - D" ,A", 



(24) 



The transition temperature and the critical field are given 
by the condition that Eq. ((20|) has a nontrivial solution 
of A" for the first time when the magnetic field and tem- 
perature decrease, respectively, where q is optimized. 



n=0 



The Abrikosov functions (/!'n(f) are expressed as 

Mi) = -l={4rMr) , (19) 



where (^o(f) denotes the solution of r]}(j)o{r) = 0. In 
the gap function (|18p . we have taken into account the 
possibility of nonzero q for the FFLO state, where q = 
{m/mxy/^\q\. The linearized gap equation (fT4|) can be 
written as 



with 



(20) 



(21) 



A. The case of s-wave pairing 

For s-wave pairing, we insert 7s (p) = 1 into Eq. (|2ip . 
It is easily verified by integrating over <f that E^^, — 
for n n' . Therefore, the matrix elements in the 

linearized gap equation (|24p are expressed as 



OO /'-TT 



ip Jo Jo sinh(2^rp/%) 



(25) 



B. The case of ci3,2_y2-wave pairing 

We consider c!a.2_j,2-wave pairing interaction given by 
Eq. ini) with 7d,2_^2 (p) = Ciivl-pl), where {px.Py) = 



5 



[PxTPy]/ \ Px '^Py- The function 7d^2_ 2 (p) rewritten 



7d.2_„2 (P) = CHzT-^ 



rrixPi + ruypl 
2mx 



nix + my tan 9 cos^ ip 



(26) 



where we have defined the polar coordinates {9, (p) tak- 
ing the px axis as the polar axis, so that 9 is measured 
from the px axis, and (p is measured from the py axis. 
Therefore, in the linearized gap equation (j24p . the ma- 
trix elements D^^r^^ are written as 



+ X! -^"n' '''''' (^n,«'+2i + '^n+2i.n') 



where -D„^^/ is defined by 



dp 



7O) (i^ 



X / singdg " — — 

Jo smh(27rTp/t;F) 



(28) 



with 



C, 



.(0) 



1 + — 



nix + my 



cjO) - 1 
c(0) + f 



x{6,,o~2ci9)+Aj[c{9)]' 



(29) 



and c(0) = 1/^1 + {my/mx)ta.Ti 9. The matrix equa- 
tion is decoupled two sets of equations for A^""^""^ 
of even and odd n's. 



C. The case of da^^-wave pairing 
Next, we consider d^^y-wave pairing interaction with 



./mxm,, tan 9 cos w 

7d.„ (P) = = \ . (30) 

mx + my tan 9 cos^ 

The matrix elements D'jj', take the form 



(31) 



where I?^^"^"*^ is defined by 



^'F Jo Jo 



with 



sinh(27rTp/%) 

(32) 



X <^ 2c 



.(0~){l-[c(0~)]2}<5,.o 



(33) 



4j 



:(0)]2J^ 



c{e) - 1 

c(^) + 1 



The order parameter components Aq""" and A'j''"'' cou- 
ple only with Aj^" and Aj^i^i.]^, respectively, where n — 
1,2,... . Therefore the matrix equation ( IMt for dx„- 



(27) wave pairing is decoupled into two sets of equations for 



An""" of even and odd n's. 



III. EXTENDED JACCARINO-PETER 
MECHANISM 

In this section, we review an extended Jaccarino-Peter 
mechanism in antiferromagnetic superconductors fzl, [§] . 
We have formulated the upper critical field equation of 
the conduction electron system in the previous section. 
Now, we shall introduce a localized spin system that is 
coupled with the conduction electron system through the 
Kondo interaction Jk. We consider a situation in which 
the antiferromagnetic long-range order has canted spin 
structures in magnetic fields. Then, ferromagnetic mo- 
ments are induced in the localized spins and create ex- 
change fields on the conduction electrons. The exchange 
fields partly or completely cancel the magnetic fields in 
the Zeeman energy terms of the conduction electrons 
when Jk > 0. The exchange fields are not real mag- 
netic fields in the sense that they only modify the Zee- 
man energy terms but not the vector potential. On the 
other hand, the ferromagnetic moments also create real 
internal magnetic fields. However, such internal mag- 
netic fields are weak when the localized spins are spatially 
separated from the conduction electrons. Therefore, we 
neglect them in this study. 

We examine the case in which the antiferromagnetic 
transition occurs at a temperature much higher than 
the superconducting transition temperature and thus 
the magnetic ordered state is rigid. Therefore, we ne- 
glect modification of the magnetic structure by occur- 
rence of superconductivity. Strictly speaking, the lo- 
calized spin state is modified by the mobile electrons, 
and the state of the total system should be determined 
self-consistently 40]. However, we consider a model in 
which such modification has already been included and 
the canted spin structure occurs as a result. 
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We modify the extended Jaccarino-Peter mecha- 
nism for the system with anisotropic Kondo interac- 
tions, introducing Kondo coupUng constants between 
the locahzed spins and the conduction electrons of the 
component in spin space. When the magnetic field 
is oriented to the —x^ direction, where /i — 1,2,3, we 
replace Hm in Eq. ^ with the effective Hamiltonian 



20 



(34) 



where denotes the effective Zeeman energy. There- 
fore, in the critical field equations derived in the previous 
section, the Zeeman energy h is replaced with h^. For 

H < Haf, is given by 



^ 10 
03 
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h=\l- 



^AF 



(35) 



with J^p = gczJ/{gsZK). In Eq. ([351), z, zk, J, and gs de- 
note number of antiferromagnetic bonds for a given site, 
number of lattice sites that participate in the Kondo in- 
teraction with the conduction electrons at r, coupling 
constant of the exchange interaction of the localized 
spins, and the g factor of the localized spins, respectively. 
Here, H^f denotes the critical field of the antiferromag- 
netic phase which is given by Haf = '^zJS / gsfJ.B, where 
S denotes the magnitude of the localized spins. There- 
fore, when < < 2J^p, the Pauli paramagnetic pair- 
breaking effect is reduced in antiferromagnetic supercon- 
ductors. 

For H > Haf, since the spins are aligned uniformly, 
the effective Zeeman energy is given by 



(36) 



Because the second term is a negative constant when 
> 0, the present mechanism is reduced to the conven- 
tional Jaccarino-Peter mechanism. The Zeeman energy 
is completely compensated at H = z^Jji^S /\fj,c\ = ^^ont- 
It is noted that FISC occurs only when H^^^^/ Haf — 



FIG. 1: Theoretical phase diagram for s-wave pairing in the 
magnetic fields parallel to the c axis. Here, SC stands for 
superconductivity. The solid and dotted curves show bound- 
aries of superconductivity in the presence and the absence of 
the FFLO state, respectively. The closed circles show the ex- 
perimental data of the superconducting transition points for 
H II c in K-(BETS)2FeBr4 by Konoike et al. (Ref. [i|). The 
dot-dashed line shows the critical field of the antiferromag- 
netic phase. 
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IV. APPLICATION TO AN ORGANIC 
SUPERCONDUCTOR 

In this section, we apply the present theory to the anti- 
ferromagnetic organic superconductor K-(BETS)2FeBr4 
taking the x and y axes along the crystallographic c 
and a axes, respectively. Therefore, we write the ef- 
fective masses as rrix = tIc, my — ma, and to^ = m?,, 
and the anisotropic Kondo interactions as = and 
= Jk. In K-(BETS)2FeBr4, the canted spin struc- 
ture was observed when H || c, and we can apply the 
extended Jaccarino-Peter theory. However, for the mag- 
netic field parallel to the magnetic easy axis, i.e., H || a, 
the canted spin structure has not been observed and the 



FIG. 2: Theoretical phase diagram for s-wave pairing in the 
magnetic fields parallel to the a axis. The notations are the 
same as those in Fig. [T] Since the canted spin structure does 
not occur for H || a at low field, the upper critical field of 
low-field superconductivity has not been calculated. 



details of the metamagnetic transition have not been re- 
vealed. Therefore, we do not apply the present theory to 
the low- field region for H || a. 

It is convenient to introduce a constant 

\mbj \2ttTPJ c 
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with which and Ky are expressed as 
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r 
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We also define the strength ratio of the paramagnetic and 
orbital pair-breaking effects 



h/2TTT, 



(0) 



(39) 



We explain the procedure to analyze the experimen- 
tal phase diagram with the present theory. Our critical 
field equations contain five microscopic parameters: Om, 



K' 



and J^. 



First, we fit the experimental 



data of the critical field of the FISC for H 



with the 



parameters a„i, rric/ma, r„i, and J^. Secondly, using the 
parameter values of a„i, rric/ma and r,„ determined in 
the first step, we fit the data of the FISC for H || a with 
a single parameter J^. It is noted that all five parame- 
ters are determined only from the curves of the FISC for 
H II c and H || a. Crudely speaking, the curve of the 
FISC can be characterized by three real numbers, that 
is, the magnetic field where the transition temperature 
is maximum, the maximum transition temperature, and 
the width of the magnetic fields at T = 0. Therefore, we 
obtain six real numbers from the two experimental data 
of the FISC, which sufficiently determine all five micro- 
scopic parameters. Lastly, using those five parameters 
determined in the first and second steps, we calculate 
the upper critical field of low-field superconductivity for 
H II c, where the canted spin structure coexists. In this 
last step, we do not use any additional fitting parameter. 
We carry out this procedure for s-wave, c?3.2„y2-wave, 

and dxy-vfave pairings. In Figs.[l]-l7l we set Tc^^ = 1.4 K 
and Haf — 4.6 T, which have been observed in the com- 
pounds K-(BETS)2FeBr4. For all pairing symmetries ex- 
amined, the theoretical curves agree well with the exper- 
imental data, as shown in Figs. [THHl In particular, the 
agreements at low field for H || c are quantitative, be- 
cause all parameters are determined only from the data 
of the FISC at high fields. In what follows, we explain 
the differences depending on the pairing symmetries. 

The results for s-wave pairing are shown in Figs. [T] 
and m We obtain the effective-mass ratio rric/ma = 
0.785, which is reasonable in comparison with the value 
estimated from the Shubnikov-de Haas (SdH) oscilla- 
tions [U |43|. It is confirmed that the FFLO state is 
not realized in the FISC. In contrast, it is shown that 
low-field superconductivity exhibits the FFLO state at 
low temperatures. 

The results for dx2-y2-waye pairing are shown in 
Figs. [3] and m We obtain rric/ma — 0.468, which is too 
small in comparison with the value estimated from the 
SdH oscillations [4l|,|4l|. It is confirmed that the FFLO 



H//c (a',5_,2-wave) 



•/k =3.oiy'Ap 

a„ = 0.0611 [T 
r„ = 0.678 
= 0.468 m„ 




T[K] 

FIG. 3: Theoretical phase diagram for dj.2_y2-wa.ve in the 
magnetic fields parallel to the c axis. The notation is the 
same as those in Fig. [T] The inset shows the theoretical phase 
diagram below 0.8 K and between 2 and 2.6 T. 
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FIG. 4: Theoretical phase diagram for d^2_y2-wave in the 
magnetic fields parallel to the a axis. The notation is the 
same as those in Fig. [T] 



state is not realized in the FISC. There are bends in the 
theoretical curves of the low fields, where the Landau 
level indexes n's of the optimum solution alternate. The 
solution near T = is expressed by a linear combination 
of the Abrikosov functions with odd n's. In Fig. O it is 
found that there are two tricritical points in the phase di- 
agram for H II c. One at T ~ 0.67 K is the critical point 
of the normal state, the BCS state with even n's, and the 
FFLO state with even n's, and another at T ~ 0.32 K is 
that of the normal state, the FFLO states with even n's 
and odd n's. 

The results for dxy-'wave pairing are shown in Figs. [5] 
and[ni We obtain rric/ma = 0.865, which is reasonable in 
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T[K] 

FIG. 5: Theoretical phase diagram for dxy-wave in the mag- 
netic fields parallel to the c axis. The notation is the same as 
those in Fig. [T] The inset shows the theoretical phase diagram 
below 0.8 K and between 2 and 2.6 T. 
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FIG. 6: Theoretical phase diagram for d^y-wme in the mag- 
netic fields parallel to the a axis. The notation is the same as 
those in Fig. [1] 




FIG. 7: Temperature dependence of the magnitudes of the 
total pair momentum along the upper critical field curves for 
(a) s-wave, (b) cij.2_j^2-wave, and (c) dxy-wa,ve pairings in the 
magnetic fields parallel to the c axis. The sets of parameters 
for (a), (b), and (c) are the same as in Figs. [1] |3l and [5] 
respectively. 



Landau-level indexes n of the optimum solution alter- 
nate. For dxy-wave pairing, it is found that the q van- 



ishes for T < 0.32T, 
n is unstable. 



(0) 



where the FFLO state with odd 



comparison with the value estimated from the SdH oscil- 
lations [4l|,|42|. In this case, the area of the FFLO state 
is very limited. It is found that the FFLO state is not 
realized in the FISC and low-field superconductivity near 
T = 0. There is a bend in the theoretical curve, where 
the parity of the Landau- level indexes n's alternates. The 
solution of the odd n's gives the highest critical field near 
T = 0. 

In Fig. [71 we show temperature dependences of g = 
(TO/mc)^^^|q| along the upper critical field curves of low- 
field superconductivity for s-wave, d^^-y^-waye, and dr^y- 
wave pairings when H || c. For dx2_y2-wsLve pairing, it 
is found that the q jumps at a temperature where the 



In a brief summary of this section, experimental results 
could be reproduced for all s-wave, dx2-y2-wa,Ye, and d^y- 
wave parings. For s-wave pairing, the upper critical fields 
of low-field superconductivity are smaller than the ex- 
perimental data, but they are improved by taking into 
account the FFLO state. The in-plane mass anisotropics 
are estimated as rric/ma = 0.785, 0.468, and 0.865 for s- 
wave, c?^2_j,2-wave, and c?a:i;-wave pairings, respectively. 
The results for s-wave and c?a;j,-wave pairings agree with 
the value obtained from the SdH oscillations [4l|, [i^ . 
whereas the result for dx2-y2-'wave pairing largely devi- 
ates from it. We found that the FFLO state does not 
occur in the FISC, while it occurs in low-field supercon- 
ductivity. 
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FIG. 8: The mass anisotropy dependence of the pure orbital 
hmits for H || c. The open circles, open diamonds, and open 
triangles show the results for s-wave, d^2 _y2-wetve, and dxy- 
wave pairings, respectively. The solid, dotted, and dashed 
curves are guides to the eye. Here, Hc20 denotes the orbital 
limit Hc20 of the isotropic case (rric = ma). 



V. SUMMARY AND DISCUSSION 

We have studied the upper critical fields of antiferro- 
magnetic s-wave, (i^2_j^2-wave, and d^y-wave supercon- 
ductors with Fermi surface anisotropy. We have de- 
rived the linearized gap equations when the magnetic 
field is applied parallel to layers, taking into account the 
effective-mass anisotropy, the FFLO state, and the ex- 
tended Jaccarino-Peter mechanism. We have applied the 
theory to the organic superconductor K-(BETS)2FeBr4, 
and obtained good agreement between the theoretical 
and experimental results. 

While the agreements have been obtained for all cases 
of s-wave, dj.2_j^2-wave, and d^y-wave pairings, the re- 
sultant in-plane mass ratios mc/ma are different from 
one another. The mass ratio for dj.2_j,2-wave pairing 
seems too small, although those for s-wave and dxy-'wave 
pairings are close to the value estimated from the SdH 
oscillations [4l|, The mass anisotropy dependence 

of the pure orbital limits of the upper critical fields is 
shown in Fig. [S] The upper critical fields for s-wave, 
(i2.2_j^2-wave, and d^jy-wave pairings become larger when 
the ratio rric/ma decreases. It is found that the upper 
critical field is the most influenced for dxy-yfSive pairing 
by a change in the ratio rric/rna, while it is less influ- 
enced for dx2 _y2-wccve pairing. This is why the resultant 
ratio rric/ma to fit the experimental data becomes much 
smaller for dj.2_j^2-wave pairing than for da;j,-wave pair- 
ing. 

This result is physically interpreted as follows. For ex- 
ample, we consider the situation that mdrna < 1. The 
origin of the orbital pair-breaking effect is the Lorentz 



force that acts on electrons. When the magnetic fields 
are applied along the c axis, the Lorentz force is the 
strongest for electrons near the Pa axis, where the Fermi 
velocity is the largest, in the two-dimensional momentum 
space. When the mass ratio mdma decreases with the 
fixed electron number, i.e., when the effective mass m,a 
increases, the Fermi velocity in the a direction decreases. 
Therefore, the Lorentz force acting on the electrons on 
the Fermi surface near the Pa axis becomes stronger and 
superconductivity is more suppressed when mdma de- 
creases. The order parameter of d^y-^ave pairing has 
nodes on the Pa axis, while that of (ia;2_j^2-wave pairing 
has the largest amplitude there. Therefore, dx2-y2-weLye 
pairing is less affected than dxy-v/ave pairing by a change 
in the ratio md ma- 
lt has been found from Fig. [5] that the FFLO state 
does not occur for dxy-wave pairing in low-field super- 
conductivity near T = 0. This result is physically ex- 
plained as follows. Since the effective density of states 
for a- wave pairing is proportional to [7q(p)]^ in the gap 
equation, areas of the Fermi surface near the peak of 
[7c((p)]^ have a greater effect on the nesting condition 
for the FFLO state, but those near the nodes of [7q(p)]^ 
have less effect \4^. Hence, in the absence of the or- 
bital pair-breaking effect, the vector q tends to points in 
the areas of the Fermi surfaces where [7q;(p)]^ exhibits 
a peak. In this study, however, we adopt the model in 
which q 11 H is assumed due to the orbital pair-breaking 
effect. Therefore, for d^y-'wave pairing, the Fermi sur- 
face nesting is less effective when H || a; or y, because 
the order parameter has nodes on the Px and Py axes. 
In contrast, for dx2-y2-wave pairing, because the order 
parameter has peaks on the Px and py axes, the FFLO 
state occurs in low- field superconductivity, when H || x 
or y. 

In the application to K-(BETS)2FeBr4, it is found that 
the FFLO state is not present in the FISC in contrast 
to A-(BETS)2FeCl4, while it is present in low-field su- 
perconductivity. In A-(BETS)2FeCl4, the phase bound- 
ary of the FISC can be well reproduced if the possibil- 
ity of the FFLO state is taken into account For 
the FFLO state to occur, the orbital pair-breaking effect 
needs to be sufficiently suppressed. In A-(BETS)2FeCl4, 
this condition seems to be satisfied even at high field 
due to low dimensionality. In contrast, it was suggested 
that K-(BETS)2FeBr4 has a larger interlayer electron 
hopping energy than A-(BETS)2FeCl4 0]. Therefore, in 
K-(BETS)2FeBr4, the orbital effect is substantial at high 
fields, as verified by comparison of the zero-field transi- 
tion temperature and the maximum transition tempera- 
ture of the FISC. 

The upper critical fields of low-field superconductivity 
agree with the experimental data, as shown in Figs. [1] [31 
and [51 Further, for s-wave pairing, the agreement is im- 
proved if the FFLO state is taken into account. However, 
they are slightly smaller than the experimental data. 
This discrepancy may be removed by taking into account 
a mixing effect in the presence of the weak triplet pairing 
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FIG. 9: Theoretical predictions of the phase boundary of su- 
perconductivity for various values of J^. The other parame- 
ters are common to Fig. [T] The thin dotted and thin two-dot- 
dashed curves show the Pauli paramagnetic and orbital limits, 
respectively. In this figure, the FFLO state is taken into ac- 
count, but modifies the curves slightly at low temperatures. 
At the upper critical fields of the FISC for = 2.5Jaf, the 
FFLO state does not occur. 



interaction, as in the compounds A-(BETS)2FeCl4 [2^ . 
A slight discrepancy also occurs in the lower critical fields 
of the FISC, where the FFLO state does not occur as ex- 
plained above. If the mixing effect occurs, this discrep- 
ancy may also be removed. 

For (ia;2_y2-wave and dxy-wetve pairings in low-field su- 
perconductivity, we found an internal transition between 
the vortex states with different Landau level indexes 
n, which is analogous to internal transitions in a two- 
dimensional system in tilted magnetic fields ^36^. Both 
the FFLO state and the vortex state with higher n's orig- 
inate from the Pauli paramagnetic effect as discussed in 
Ref. [1^, because in the absence of Zeeman energy only 
the vortex state with lower n's occurs. Due to the transi- 
tions to the FFLO state and the vortex states with higher 
n's, the upper critical fields of low field superconductivity 
exhibit downward convex curves for all pairing symme- 
tries examined, which agrees with the experimental data. 

We have determined the microscopic parameters only 
from the curves of the FISC, and calculated the upper 
critical field of low-field superconductivity with those pa- 
rameter values without any additional fitting parame- 



ters. From the accordance between the theoretical and 
experimental results on low-field superconductivity, the 
extended Jaccarino-Peter mechanism 0, ll] seems to be 
realized in K-(BETS)2FeBr4. Applying the extended 
Jaccarino-Peter mechanism to the same models with dif- 
ferent parameter values, we obtain some unusual phase 
diagrams as shown in Fig. [HI When = l.SJ^p, the 
transition temperature curve has double peaks. When 
Jk — J'kfi superconductivity occupies large area in the 
phase diagram. Similar phase diagrams have been ob- 
tained in two- and three-dimensional systems [1, [2^ . In 
this study, we confirmed that the phase diagrams as ob- 
tained above can be realized in a more realistic model 
of organic superconductors. If we control the value of 
Jk by material design, such critical field curves would be 
observed. 

In the present theory, the effective mass m;, is included 
in the parameter a^. Since we have obtained the value 
of arm we can estimate a value of m?,, if we have the value 
ofpF appropriate for K-(BETS)2FeBr4. 

The upper critical field in low fields for H || a could 
not be reproduced by the present theory. It has been 
found that the magnetic easy axis of antiferromBgnetic 
long-range order in low fields is along the a axis [10| , and 
that the metamagnetic transition occurs at the antifer- 
romagnetic phase boundary close to the upper critical 
field [3|. In order to reproduce the upper critical field in 
low fields for H || a, we need to reproduce the localized 
spin structure. This remains for future study. 

In conclusion, the experimental phase diagrams of the 
antiferromagnetic superconductor K-(BETS)2FeBr4 are 
theoretically reproduced by the present theory. In partic- 
ular, the low-field phase for H || c is well reproduced by 
the model parameters that are determined from the criti- 
cal fields of the FISC. Therefore, the extended Jaccarino- 
Peter mechanism seems to be realized in the present com- 
pound at low fields. Due to this mechanism, some un- 
usual phase diagrams may occur in compounds with ap- 
propriate energy parameters, in the future. 
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